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Abstract 



The chemical evolution of nascent quark matter core in a newborn 
compact neutron star is studied in presence of a strong magnetic field. 
The effective rate of strange quark production in degenerate quark 
matter core in presence of strong magnetic fields is obtained. The in- 
vestigations show that in presence of strong magnetic fields a quark 
matter core becomes energetically unstable and hence a deconfinement 
transition to quark matter at the centre of a compact neutron star 
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under such circumstances is not possible. The critical strength of mag- 
netic field at the central core to make the system energetically unstable 
with respect to dense nuclear matter is found to be ~ 4.4 x 10 13 G. This 
is the typical strength at which the Landau levels for electrons are pop- 
ulated. The other possible phase transitions at such high density and 
ultra strong magnetic field environment are discussed. 
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1. INTRODUCTION 



More than a decade ago, it was suggested by Witten [p]] that a flavor symmetric 
mixture of u, d and s-quarks known as strange quark matter (SQM) would be the 
most stable configuration of hadronic matter. An extensive theoretical investigations 
on the stability of both bulk as well as quasi-bulk properties of SQM using MIT bag 
model and also with some other confinement models have proved the reality of 
such speculation [@-@]. Usually the strange quark is attributed a minor role in the 
standard model of the Universe. However, it has also been suggested by Witten that 
strange quark matter might be an important relic from the early Universe and also 
a constituent of neutron stars. Although strange quark matter could be absolutely 
stable and thus be the true ground state of hadronic matter, long-lived or stable 
strange quarks are only found in the regions with very high density, otherwise, the 
relatively large s-quark mass and the weak mixing to the down quarks lead to rapid 
decay. In a quark matter system, either at the core of a compact neutron star or 
in the early Universe-microsecond after big bang, s-quarks are produced through 
a number of weak processes. It is generally believed that microsecond after big 
bang, the Universe was completely filled with a hot and dense soup of quark matter. 
The survival of strange quark nuggets relic till today from the early Universe is a 
controversial issue ||. It is also expected that a deconfinement transition to bulk 
quark matter may cur at the core of a neutron star if the density exceeds a few times 
normal nuclear matter density. In such a cold quark matter system, s-quarks are 
produced through weak processes. Although the existence of strange quark nuggets 
relic produced as a result of first order QCD phase transition in the early Universe 
is highly uncertain, the presence of exotic quark matter core in a compact neutron 
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star can not be ruled out. A few years ago, Madsen and Hiselberg et. al. have 
investigated with detail numerical calculations the rate of s-quark production in 
dense quark matter in the context of survival of strange quark nuggets produced 
at the time of primordial QCD phase transition |f|[7j. Recently, Glendenning and 
his group have reported a lot of exciting results on the existence of strange quark 
matter core in a compact neutron star and its effect on various gross-properties 
of hybrid stars ||. It has also been suggested that a transition to quark matter 
core may occur during slowing down of a old neutron star ||. A first order phase 
transition to quark matter at the core of a neutron star is also possible by seeding 



strange droplets from outside and by neutrino sparking ||10|| . However, in none of 
these papers, the effect of strong magnetic field of the neutron star (particularly 
in the case of newborn neutron star) on such exotic matter have been considered. 
Large magnetic fields are known to be present in neutron stars [|TT|-|i"5||. The direct 
evidence from pulsar observations show that the strength of this magnetic field at 
the neutron star surface is 10 12 — 10 14 G. Then it can very easily be shown by the 
scalar virial theorem that magnetic field strength at the core region of a neutron star 
may reach up to 10 18 — 10 19 G [14]. On the other hand, from the stability criteria 
of a compact magnetized neutron star it can very easily be shown following the 



famous work by Chandrasekhar and Fermi [[L5[] in the context of magnetized white 
dwarfs, that the upper limit for the strength of magnetic field at the core region is 
~ 10 19 G. The presence of such strong magnetic field inside newborn neutron stars 
significantly modifies the order of qaurk-hadron phase transition at the core. It 
has been shown that a first order quark-hadron phase transition initiated by the 
nucleation (by thermal or quantum fluctuation) of quark droplets in a system of 
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meta-stable neutron matter in presence of a strong magnetic field (> B$^) is 
absolutely forbidden P,|T6||. However, some higher order transitions, e.g. metal- 
insulator type of second order phase transition is allowed if the strength of magnetic 
field is < 10 20 G. Which is of course too high to achieve at the core of a neutron 
star. Even if the transition is of second order in nature, the gross properties of 
dense quark matter are considerably modified with respect to the field free case. 
The presence of strong magnetic field also affects the weak and electro-magnetic 
processes which may occur in the degenerate quark matter core |JL7|| . The presence of 
strong magnetic field shifts the /3-equilibrium point of the quark matter system |18| 



It also affects the abundances of various components in /3-equilibrium. The photon 
emission by bremsstrahlung, photon splitting and neutrino emission processes are 



significantly affected by the presence of strong magnetic field [T7JTS]. It is also 
believed that pulsar kicks are generated by the asymmetric neutrino emission in 
presence of strong magnetic fields |f20|-p3|]. In a dense quark matter, since s-quarks 
are produced or annihilated through weak processes, the presence of strong magnetic 
field should also alter the production and decay rates of s-quarks. To the best of 
our knowledge, theoretical investigation of this particular problem has not been 
reported earlier. In the previous work by Madsen and Hiselberg et. al. ||[7| the rate 
of s-quark production in degenerate cold quark matter (T < lOMeV) was studied 
for the field free case in the context of stability of strange nuggets as the relic of 
big bang deconfinement transition. In the present paper we shall study the effect of 
strong magnetic field on the net rates of u, d and s-quarks production and there by 
study the chemical evolution of nascent degenerate quark matter core in a compact 
neutron star. We shall also consider the effect of neutrino trapping in degenerate 
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quark matter on the chemical evolution in presence of a strong magnetic field. The 
main finding of this paper is that the chemical evolution of nascent quark matter core 
in presence of a strong magnetic field leads to an unphysical scenario- the final state 
in /^-equilibrium is energetically unstable with respect to the normal hadronic matter 
of neutron star. Although a flavor symmetric mixture of u, d and s-quarks at zero 
temperature is energetically much more stable than the corresponding field free case, 
the chemical evolution of a bulk system consisting of u, d and s-quarks (not a flavor 
symmetric mixture) shows that the quark matter core of a newborn neutron star 
may become energetically unstable if the magnetic field strength exceeds a typical 
value ~ 4.4 x 10 13 G, which is the strength of magnetic field at which the Landau 
levels for electrons are populated and quantum mechanical effect of strong magnetic 
field on electrons become important. Since the formation of a flavor symmetric quark 
matter at the central region of a compact neutron star is possible only through the 
chemical evolution of ordinary quark matter or u, d and s flavor asymmetric quark 
matter, it would really be interesting to investigate the time evolution of various 
species in the system both for field free case and in presence of a strong magnetic 
field. 

For the sake of simplicity we neglect quark-quark interaction in the degenerate 
quark matter. The temperature of the system is taken to be lMeV (~ 10 9 K). 
The increase or decrease in temperature of the system causes an enhancement or 
suppression of effective production rate of s-quark respectively. However, it only 
makes a quantitative change in our result. The qualitative nature remains unaltered. 
It is further assumed that the external magnetic field B m is constant throughout the 
matter. The convenient choice of gauge is A = and A = (0, xB rn , 0) corresponding 
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to the constant magnetic field B m along z-axis. If the strength of magnetic field B m 
is greater than the critical value Bffi® = mf/qi for the i-ih charged particle, where 
rrii and qi are respectively the mass and charge of the i-th species, the corresponding 
single particle energy is given by £j = (p 2 z + m\ + 2q i uB m ) 1 ^ 2 , otherwise the usual 
expression for field free (B m = 0) case is used, where v is the Landau quantum 
number. The critical strength of magnetic field as mentioned above can be obtained 
by equating the cyclotron quantum with the rest mass energy of the i-th charged 
particle. If the magnetic field strength exceeds this typical value, the quantum 
mechanical effect of strong magnetic field on the i-th charged component becomes 
important. This is known as Landau diamagnetism. For electrons of mass 0.5MeV, 
this typical strength is B^^ ~ 4.4 x 10 13 G, whereas for u and rf-quarks of current 
mass 5MeV, it is ~ 10 2 Bffl e \ On the other hand for s-quark of current mass 
150MeV, this typical value is ~ 10 20 G. Which is too high to achieve even at the core 
of a newborn neutron star. Therefore in the case of s-quarks, quantum mechanical 
effect of strong magnetic field is insignificant and results for the field free case are 
used. 

The paper is organized in the following manner, in section 2 we have given the 
basic formalism to study the chemical evolution of compact and degenerate nascent 
quark matter core with non-degenerate neutrinos (i.e. neutrinos are assumed to 
leave the system freely immediately after their formation), both in absence and 
presence of strong magnetic fields. In the case of B m ^ 0, we have considered two 
possible scenarios- (i) magnetic field strength is low enough so that only electrons 
are affected and (ii) the field strength is relatively stronger so that all the charged 
particle components except s-quarks are affected. In section 3 we have discussed with 
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detail numerical calculations, the effect of neutrino trapping on chemical evolution 
of the system for both B m = and B m ^ cases. Concluding remarks are presented 
in section 4. 

2. BASIC FORMALISM: WITH NON-DEGENERATE NEUTRINOS 

To study the chemical evolution of degenerate quark matter at the core of a 
newborn neutron star we consider the following weak interaction processes which 
are taking place at the quark matter core 

d — > u + e~ + v e (1) 

u + — > d + u e (2) 

s — > u + e~ + V e (3) 

u + e~ — > s + v e (4) 

■u + d # m + s (5) 

We further assume that the system is charge neutral and produced neutrinos/anti- 
neutrinos leave the system freely (i.e. they are non-degenerate). The approach 
to chemical equilibrium of the system is governed by the following sets of kinetic 
equations 

^_L| ri -r 2 + r 3 -r 4 ] (6) 
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dY d 



1 



[-r\ + r 2 - r 5 + r 6 ] 



(7) 



dt 



Baryon number conservation and charge neutrality conditions relate strange quark 
and electron abundances with that of u and rf-quark and are given by 



where Yj = rii/ns is the abundance for the i-th component (i = u, d, s, e and v e (for 
the degenerate neutrino case only)), % is the baryon number density, is the rate 
of i-th weak interaction process (i = 1, 5) and T 6 is the rate of reverse process as 
indicated by eqn.(5). 

To study chemical evolution of the system we solve the differential eqns. (6) and 
(7) numerically along with two constraint equations, eqns. (8) and (9) and use the 
rates r\ to T 6 , appear on the right hand side of eqns. (6) and (7). Since these 
differential equations are of first order in nature, we need at least one initial condition 
for both Y u and Y d in order to get the complete solutions. For a neutron star of 
mass ~ 1.4M the baryon number density % at the centre is of the order of 3 — 4no, 
temperature ~ 10 9 K and proton content is about 4%. Then the initial conditions 
are given by Y u (t = 0) = 1.04, Y d (t = 0) = 1.96. As a consequence of baryon number 
and charge conservation we have Y s (t = 0) = and Y e = 0.04. 



In this subsection we shall first derive the rates (I\ to T 5 ) for the processes given 
by eqns.(l)-(5) and also the rate T 6 for the reverse process as indicated by eqn.(5), 
for the field free case. 



Y s — 3 — Y u — Y d 



(8) 



Y e = Y u - 1 



(9) 



1. Weak Reaction Rates: B m = 
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The reaction rate for the process (1) is given by 

r\ = 6 x 64G 2 cos 2 9 c jj2 jS^T^^ 1 ~ & ~ ** ~ P*) 

/( £l )(l - /(e s ))(l - /(£4))(pi.p 2 )(p 3 .p4) (10) 

where z = 1,2,3,4 stand for d, v e , u and e respectively, /(£;) = 1/(1 + 
exp((ej — fjLi)/T), the equilibrium (Fermi) distribution function for the z-th com- 
ponent and Pi, £j, /ij are respectively the 4- momentum, energy and chemical po- 
tential of the particle i, G is the weak coupling constant, 9 C is the Cabibbo angle 
(cos 2 6 C ~ 0.74) and T is the temperature of the system. The angular part of the 
integrals can be evaluated without any approximation |||7| and we get 

T 1 = — — cos 2 B c / 53 deiS(ei - e a - e 3 - e 4 ) 
71 J i=i 

/( £l )(l - /(e 3 ))(l - /(e 4 ))/„(ei,e 2 ,e3,e4) (11) 

where 

I n (ei, e 2 , S3, e 4 ) = &i 2 A;+ £2, £3, £a) + ~{k+ - k l })I 2 (ei, e 2 , e 3 , e 4 ) - 



— I 3 (£i,£2,£3,£4) (12) 
^ = e iei ±i(p?+pj) (13) 

HPS 1 ' 1 - PtT^Pk - PjMpm - PiMPij - Pkl ) 

HPtr 1 - pir 1 )^ - Pk)o(p u - Pa) 

HPlr 1 - PS' 1 )^ - P*)0fa - P»)0(Pk - Pa) (14) 
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Pij = Pi + Pj, Pij =| Pi — Pj \, Pi —\ Pi | is the magnitude of three momentum vector 
and Si = (pf + m 2 ) 1//2 (we assume h = c = k B = 1). 

In the case of a compact neutron star, the central density is very high (~ 3 — 4n 
to have a deconfmement transition) but the temperature is low enough (~ 10 9 K), as 
a consequence u, d quarks and electrons are relativistic in nature and are strongly 
degenerate. Under such an extreme condition all the integrals except one in eqn.(ll) 
can be evaluated analytically and the final expression for rate of the process (1) is 
given by 

3G 2 r°° (x — £ d ) 2 + 7r 2 

Ti(fi d , fiu, He,T) = ^ cos2 ^y o dx - - ^ _ I(fi d , Tx, /i n , fi e ) (15) 

where in the expression for / we have replaced energy and momentum of all the 
particles except the neutrino by their values on the Fermi surface, i.e. we put 
£j = Hi and pi = PFi = (f^j — rrf) 1 ^, we have also changed the integration variable 
e v to x, where neutrino energy e v = p u = xT and £ d = (/i d — /i u — /i e )/T which 
is zero in the /^-equilibrium condition. Following the same methodology, we have 
obtained the rate for the process (2), given by 

r 2 (//d, fJ-u, He, T) = cos 2 9 C T 3 f dx j- X + ^ + W - J(n d ,Tx,n u ,n e ) (16) 

4t7T JO 1 -p Cxp yX ~\~ <,d) 

where 

J(e 1 ,e 2 ,e 3 ,e A ) = k 12 Kl 4 Ii(ei,e 2 ,e 3 ,e 4 ) 

-\( k + + k+)I 2 (ei,e 2 , e 3 , e 4 ) + ^h(ei, ^, e 3 , e 4 ) (17) 

The reaction rates for the processes (3) and (4) corresponding to s-quark can very 
easily be evaluated by replacing all the parameters for o?-quarks by the corresponding 
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s-quark values (i.e. fi d — > fM s , (,d — > £s, — ► rn s ) and cos 2 6^ by sin 2 9 C in the 
expressions for Ti and Y% respectively. 
For the weak process (5) we have 

9G 2 ej 9 „ r<jto-m.)/T 



r 5 (p d , fi u , (jL a ,T) = — — cos 6> c sin 6» C T / <ix 

Z7T° J-oo 

(x + e s ) 2 + vr 2 



J(fi d ,fi u ,fi u ,lj, 8 -Tx) (18) 



(1 + exp(ac))(l + exp(-x !E - £ a )) 
where ^ = - /i s )/T. 

The reaction rate for the reverse process as shown in eqn.(5) is related to the 
forward rate T 5 by 

r 6 = exp(-£ s )r 5 (19) 



2. Weak Reaction Rates: B m ^ 

We shall now evaluate the rates Ti to T 6 in presence of strong magnetic fields. 
We consider two possible cases: the magnetic field B m is not very high (B^^ < 
B m < B!$( U > ( V) so that the quantum mechanical effect is important for electrons 
only, and finally the magnetic field is strong enough to affect all the charged particle 
components except s- quarks (B^ u ^ < B m < B$^) . 

Let us first consider case 1. The magnetic field B^^ < B m < B$( u ' d \ To eval- 
uate the rates for the weak processes (l)-(4), we use the modified form of spinor so- 



lution of Dirac equation in presence of a strong magnetic field for electron [16,24,25 
For other charged particles, the usual plane wave Dirac spinor solutions for field free 
case are used. Since only electrons are affected, the rates for the process (5) and its 
reverse process will remain unchanged. 
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In this case the rate of the process (1) is given by 
3G 2 q e B m 4 2 ^ (pf u \ 

Tl = =-= T COS C }^ HufXe 

27r6 \PFe J 



I ( (x u + x e - ^ U + ^ — — ) f(x u )f(x e )dx u dx e (20) 

J—oo J—oo V 1 J 

where g e is the charge carried by electron, f(xi) = 1/(1 + exp(xj)), and pp e = 
(yLig — m 2 — 2vq e B m ) 1 l 2 is the Fermi momentum for electron. The sum v is over 
Landau quantum number for the electron which in principle can have all possible 
positive integer values including zero. 

Similarly the rate for the process (2) is given by 
„ 3G 2 



2tt 6 



(g eJ B m )T 4 cos 2 6 C HuHe [ — J / / dx u dx e 

u=Q \PF e ) J-oo J-oo 



x u + x e + ^ u + ^ — —^j f(x u )f(x e ) (21) 

As before the rates for the processes (3) and (4) may be obtained by simply replacing 
all the parameters for <i-quarks by the corresponding s-quark values and cos 2 8 C by 
sin 2 9 C in eqns.(20) and (21) respectively. Whereas the rates for the direct and the 
reverse processes given by eqn.(5) remain unaltered. 

Next we consider the second case possibility, B^ > B m > B$( u > d \ In this case 
modified spinor solutions for all the charged particle components except s-quark are 
used to evaluate rates of the weak processes As we have noticed that 

compact analytical expressions for the rates can only be obtained if we assume that 
zeroth Landau levels are populated for the magnetically affected charged species. 
Otherwise we have to evaluate the rates numerically right from the beginning. 

The rate for the process (1) is given by 



G 2 , „ s \2q u B m q d B r 1 1 



Ti — 7r-^{<leB m) 



PF u PF e 
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+ 00 f + OO 



/+OO r 
-oo J —oo 



x u -\- x e -\- 



Hd - Hi 



f (•£«) / (^e) dx u dx e 



(22) 



where g n and g rf are respectively the charge carried by u and rf-quarks. 
Similarly the rate for the process (2) is given by 

"2quB m qdB r , i 1 2 



T 2 = 

C+oo r+oo 



cos 2 9 r T 4 



/+oo r 
-oo J —oo 



x u -\- x e -\- 



Hu + He ~ 



/ ( x u) f ip^e) dx u dx e 



(23) 



In the present scenario, the rates for the processes (3) and (4) can not be obtained 
by simply replacing <i-quark parameters with the corresponding s-quark values and 
cos 2 9 C by sin 2 # c in eqns.(22) and (23) respectively. Since s-quarks are not affected 
by the magnetic field, we have to evaluate T 3 and T 4 with zero field plane wave 
solution for s-quark and modified Dirac spinors for all other charged species. Then 
we have the expressions for rates of the processes given by eqns.(3) and (4) for zero 
Landau quantum number 

5G 2 ■ - r+oo /-+oo 



16tt 7 



/+oo r+oo 
/ f{Xu)f{Xe) 
-oo J —oo 



x u -\- x e -\- 



Hs Hu He 
f 



dec n dec g 



(24) 



and 



5G 2 
16tt 7 



/+oo r + oo 
/ f(x u )f(x e ) 
-oo J —oo 



Hs Hv 



He 



X u ~\~ X e 

respectively. 

The rate for the process (5) is given by 



dec doc g 



(25) 



9G 2 
40tt 6 



/ d\2 2n ■ 2 a ^ HuHd 

{q e B m ) cos 0c sin 9 C \ - 



6 PF u PF d 

14 



T 1 - 



exp( 



k )/ 



f^uf^d / Jx u =-oo 
^ 2 \ / ^ 2 



dx u f (Xu) f [ X u X s + ) 



-P s Jexp( 



2q e B Tl 



PsJ 



T 

(PF U +PF d -Ps z )Ps 



5q e B n 

where d 3 p s = p 2 s dp s sin 9d9d<p. Since i? s = wp^ + m 2 s) we have 



<i 3 p s = p s E s dE s sin 9d9d(p 



f(x s )d 3 p s (26) 



(27) 



Ps* = Ps sin 9 cos 



(28) 



p Sy = p s sin 6 1 sin i 



and 



Ps z = Ps cos 9 



(29) 



(30) 



To evaluate this integral we use x s = (fx s — E s )/T as the new integration variable 
instead of E 3 . Then we have 

T-i / D \2 2 n ■ 2 n ^ f^uf^d 

r 5 = -r-^(q u B m ) cos f c sin t^ c W-- 



40tt 6 



T 1 — 



®PF u PF d 



PF u PF d 



\ Hu^d 

f(x u )f(x s )f (-x u -x s + 
( 3 , \ / 



/+oo r+oo rir r2ix 

j dx u dx s / sin#<i# / <i0 
-oo J-oo JO JO 



Hd + /is 



T 

3 



(PF U +PF d ~P Sz )Ps 



exp 



exp 



-pL 



(31) 



5q e B m Sx J 1 y 2g e .B r] 

The numerical evaluation of this mult i- dimensional integral can be made little bit 
easier (less time consuming job in computer) by assuming E 3 ~ /i s and p s ~ p Fs in 
those parts of this integral where 9, <fi and E 3 or p 3 are coupled together. The rate 
for reverse reaction is as usual given by eqn.(19). 
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To study chemical evolution of the system, we have numerically solved the sets 
of kinetic equations given by eqns.(6) and (7). As has already been mentioned, we 
have considered three possible scenarios: (i)B m = 0, (ii)B^^ < B m < B^ u ^ 
and {iii)B^ u ^ < B m < £<£)(«). The initial conditions are Y u (t = 0) = 1.96 
and Y d (t = 0) = 1.04 for the fractional abundances of u and <i-quarks. The rates 
of various processes (eqns.(l)-(5)) are also evaluated numerically and finally the 
charge neutrality and baryon number conservation conditions (eqns.(8) and (9)) are 
used to relate Y s (t) and Y e (t), the fractional abundances for s-quark and electron 
respectively with that of u and d quarks. In fig.(l) we have plotted time evolution 
of the fractional abundances Yj for i = u,d, s, and e in the field free case. This figure 
shows that ultimately fractional abundances for u, d, s-quarks and electron saturate 
to their equilibrium values. Which is the chemical equilibrium condition among 
the participants. We have further noticed, that as the current mass m s for s-quark 
approaches the current mass of u and ci-quarks, the equilibrium values for fractional 
abundances of all the three types of quarks ultimately saturates to a constant value 
1 and the electron abundance Y e m 0. We have also investigated the sensitivity 
of initial conditions. We have considered a wide variation of proton fraction- from 
almost pure neutron matter (proton fraction x ~ 0.001) to symmetric nuclear matter 
(x = 0.5). The qualitative nature of our results do not change with the change of 
initial condition. The change in initial condition only causes a change in saturation 
time scale. Now we would like to see whether such theoretical observation supports 
the speculation of Witten. We have verified whether the system in /^-equilibrium 
is energetically stable or not. To do that, we have followed the following steps, (i) 
Using these equilibrium values for various components (i = u, d, s, and e) we have 
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obtained the equation of states for the quark matter system with B m = 0, (ii) using 
the equation of state we have solved TOV equation and obtained the mass of the 
quark core, which is ~ O.4M for the central density 4n and the corresponding 
radius is R = 3.8Km. To obtain the mass of quark core in the field free case, we 
assume that the deconfinement transition is first order in nature. The corresponding 
critical density for a typical bag parameter value B 1 ^ = 160MeV is n c « 2.1n . It 
has been observed that the mass of the core increases in presence of a magnetic field 
where as the radius of the core decreases, i.e. the core becomes more compact in 



presence of strong magnetic fields [26]. To investigate the stability of such a bulk 



quark matter, we have calculated the energy per baryon using the relation 

e 1 



n B n B 



c mag. 

(32) 



where q is the internal energy density of the species i and e mag . = B^/ (8ir ) is the 
magnetic energy density, which is zero in the field free case. Using the saturation 
values for abundances for various components, we have seen that < 6qm/tib < 
M n < eNAi/ns- Which indicates that the chemical evolution of the system in field 
free case leads to a scenario which is energetically more stable than nuclear matter 
and the binding energy of the system is negative. Therefore, the observation is in 
agreement with the speculation of Witten. In fig. (2) we have plotted the variation 
of fractional abundances for various species with time for B m = 10 14 G, so that only 
electrons are affected. In the present scenario the deconfinement transition from 
hadronic matter to quark matter can not be first order. We assume the transition is 
metal-insulator type second order in nature. In this particular case, unlike a mixed 
phase of quark matter and hadronic matter in a first order quark hadron phase 
transition, we have a sharp interface between nuclear matter crust and quark matter 
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core. At the sharp boundary all the thermodynamic variables change continuously, 
i.e. P NM = Pqm, e NM = e QM , fi p = 2/j, u + fi d and fjL n = /i u + 2fi d . From these 
relations, we obtain the variation of critical density with the strength of magnetic 
field. In fig. (3) we have shown this variation. The corresponding bag parameter 
also changes with the strength of magnetic field. The dashed curve in fig. (3) shows 
the variation of B 1 ^ with the strength of magnetic field. In both the curves, up to 
B m = 10 3 B$^ only electrons are affected beyond which we have considered the 
quantum mechanical effect of strong magnetic field on all the charged species except 
s-quark. In the nuclear matter sector, we have considered a a — u exchange type of 



mean field model |24|. These curves show that both the quantities remain almost 
constant (equal to their values for field free case) up to B m = 10 17 G, beyond which 
they become extremely large. As it is obvious from fig. (2) that in this particular 
case, the net production of s-quark starts after 10 _14 sec, which is much earlier than 
the field free case, but Y s goes to zero by 10~ 10 sec and then it remains zero for 
ever. It is also obvious from the figure that in this case d quark abundance also 
goes to zero by the same time scale and in the chemical equilibrium condition it 
remains zero. The system in chemical equilibrium is therefore mainly consists of 
w-quarks and electrons. If we translate this scenario to the case of a neutron star, 
this is equivalent to a proton rich core. As a consequence direct URCA process 
will play the main role in the cooling process by neutrino emission from this region, 



provided it is energetically stable [Z7||. To study the stability of such abnormal 
system, we have calculated internal energy per unit volume. Near nuclear density 
it is ~ 4.39 x 10 14 MeV 4 , which corresponds to free energy per baryon >> M n 
and also >> energy per baryon for hadronic matter. The magnetic energy density 
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is negligibly small compared to internal energy density of the system unless the 
magnetic field is > 10 19 G. In this particular case, the overall energy of the system is 
therefore equal to the kinetic energy of the system. This indicates that the system is 
therefore energetically unstable if the magnetic field affects electrons present in the 
system. In fig. (4) we have plotted the abundances for various species against time for 
B, m = 5 x 10 16 G. In this case all the charged particles except s-quarks are affected. 
We have seen that s-quark production is absolutely forbidden in this particular 
scenario. The system behaves like a neutron rich (rich in <i-quark) core. In this case 
free energy per unit volume is ~ 4.9 x 10 12 MeV 4 near normal nuclear density, which 
again corresponds to energy per baryon >> M n and also >> energy per baryon for 
hadronic matter. Therefore, in this case also, the system is energetically unstable. 
We have also noticed that the presence of a small percentage of A-hyperons in the 
core region before phase transition does not change the conclusions. In figs. (5) and 
(6) we have plotted the abundances for the species with time for the presence of 
~ 25% A-hyperon before phase transition to quark matter for B m = and B m = 
5 x 10 16 G respectively. We have noticed that for the field free case, the time evolution 
leads to a stable configuration, where as for the presence of strong magnetic field 
the system becomes energetically unstable in the chemical equilibrium condition. 
Therefore the presence of A-hyperons in hadronic matter (which is most likely if 
the density of hadronic matter exceeds a few times normal nuclear density) before 
deconfinement transition does not change our conclusion. Next we shall consider 
the effect of neutrino trapping on the stability of quark matter core in presence of 
strong magnetic fields in the next section. 
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3. BASIC FORMALISM: WITH DEGENERATE NEUTRINOS 



In this section we shall develop a formalism to investigate the chemical evolution 
of the system considering the effect of neutrino degeneracy on the weak reaction rates 
for the field free case and also in presence of a strong magnetic field. Assuming 
the formation of a neutron star immediately after supernova explosion, then at 
the beginning the matter within the remnant, particularly at the central region is 
extremely lepton rich (neutrino trapped). Such an object is called proto-neutron 
star. Normally it takes a few second to convert a proto-neutron star to a neutron 
star by the emission of trapped neutrinos. If the density of such matter exceeds 
3 — 4no, a deconfinement transition to neutrino trapped quark matter can occur at 
the core region. To study the chemical evolution of this neutrino trapped quark 
matter, we have to evaluate the rates for the processes 

d + v e <-> u + e~ (33) 

s + u e <-> u + e~ (34) 

and the direct and reverse reactions given by 

u + d <-> u + s (35) 

In the case of trapped neutrinos in a degenerate quark matter we have to solve the 
differential equations 

^ = --V 1 + r' 2 ] (36) 

dt riB 

and 
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^ = iVi-r' 3 ] (37) 

at ub 

We also have to consider the total lepton number conservation, given by 

Y x = Y e + Y v ~ 0.4 (38) 

The magic number 0.4 is conventionally used in the numerical computation of super- 
nova explosion. The baryon number conservation and charge neutrality conditions 
are also used to obtain fractional abundances of all the species present in the system. 
In eqns.(36) and (37), T-'s are the rates of the processes (33)-(35). 



1. Rate of Weak Processes: B m = 



In the case of degenerate neutrinos, the neutrino chemical potential ji u ^ and 
has to be evaluated from eqn.(38). The rates for the processes given by eqns.(33) 
and (34) are given by 

r; = G2c ^ 9c T 3 e d (^ 2 + a,, *0 (39) 

and 

3G 2 . ,„ f^s-m s )/T 

ax 



r' 2 = ^sin^ c T 3 (l-exp(-0) 



(l + e*)(l + e-(x + ti s )) 
Finally the effective rate for the process (35) is given by 

I* = ^G 2 cos 2 6 C sin 2 9 c T\l - exp(-O) 

rO*.-m.)/T ( x + £)2 + n 2 



J{n u -Tx,n e ,n s ,n u ) (40) 



/(Hs-m s )/i i x _|_ f y _|_ Tj-" 

^ (l + e^(l + e-(^)) J ^ ^ - Tx ' ^ **) ( 41 ) 
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2. Rate of Weak Processes: B m / 



Now we shall consider the weak process in the degenerate quark matter in pres- 
ence of a strong magnetic field and also consider the effect of neutrino neutrino 
trapping. As before, we first assume that B^^ < B m < B^ u ' d \ so that only 
electrons are affected quantum mechanically by the magnetic field. In this case the 
effective rate for the process (33) is given by 

1 



G 2 cos 2 B c ~ (p F \ r°° yoo 

T 1 = (Tq e B m )}^ fi u fj, e / dx u / dE 2 



[exp(-^/T - x u ) + 1] 
E'J(x u )[l + exp(-(£ 2 - n v )/T)\ (42) 

where E' 2 = E 2 - (fx u + (J, e ~ Hd)- 

To obtain rate for the process (34), we have to replace all the parameters of 
(i-quark by the corresponding s-quark parameters and cos 2 9 C has to be replaced by 
sin 2 9 C . In this case E' 2 = E 2 — (fi u + A* e — aO- The rate for the process (35) is given 
by eqn.(41), which is the expression for field free case. 

When the magnetic field strength B m > B^ u ^ but less than B^ s \ all the 
charged species, except s-quarks are affected by the strong magnetic field. As has 
been mentioned before, compact analytical expressions for rates can only be obtained 
if we assume that only the zeroth Landau levels are populated for all the magnetically 
affected charged particle components. In this case the rate for the process (33) is 
given by 

3G 2 / q u B m qdB m \ ^ A*uA*e 



■p/ 3G / q u B m q d B m \ 2 fJ> u fj, e 

r i = T^VeBm — cos 6 C T 

2vr 4 V 17 J p Fu p Fe 

Jx e =-oo Je 2 =o Je=o J<t>=o \ fidE 2 J 

[1 + exp(-(£ 2 - ^)/T)]Eij— — - ^^f( Xe )dx e dE2 sm9d0d<f) (43) 

[1 + exp[ x e cj 2 I 1 )\ 

22 



As before the rate for the process (34) can not be obtained simply by replacing 
ci-quark parameters by the corresponding s-quark values and cos 2 9 C by sin 2 9 C . In 
this case the rate is given by 

r' 2 = ^(«A.ftA.) (it) (it) 1*1 sin* C (44) 

I = J J (x u + x e - — — y — — ) f(x u )f(x e )dx u dx e (45) 

The rate for the process (35) is given by the combination of eqns.(31) and (19). 

In the degenerate neutrino case, we have plotted the time evolution of fractional 
abundances for various species in the filed free case in fig. (7). In this case the 
system is marginally stable, internal energy density is ~ 6.25 x 10 9 MeV 4 . Since 
it takes a few second to convert such a system to conventional system of SQM 
with non-degenerate neutrinos by neutrino emission, we may conclude that it is 
possible to have an energetically stable quark matter core from the deconfinement 
transition at the central region of a compact proto-neutron star. In the case when 
< B m < B$( u ' d \ i.e., only electrons are affected, we have seen that al- 
most all the neutrinos are absorbed in the system momentarily and to keep lepton 
number constant, electron density increases by several orders of magnitude. As a 
consequence of charge neutrality condition, the w-quark density also increases by the 
same orders of magnitude. But the conservation of baryon number of the system 
causes decrease in rf-quark density almost by the same orders of magnitude. The 
situation is happened to be the most unphysical one. In this case the binding energy 
calculation shows that the system is extremely unstable. In this case internal energy 
density is ~ 10 18 MeV 4 near normal nuclear density. In the case of strong magnetic 
field (B^ u ^ < B m < B$^), when all the charged components except s-quark are 
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affected, we have shown the time dependence of fractional abundances in fig. (8). In 
this case also s-quarks are never produced and the conclusion does not change if 
there are a few A-hyperons at the central region of proto-neutron star before phase 
transition. We have also seen that electron density goes up by several orders of 
magnitude, as a consequence -u-quark density also increases and the system behaves 
like a neutrino trapped proton core (if we translate it to proto-neutron star without 
quark core). The scenario is also unphysical in this case. The internal energy of the 
system is 4.376 x 10 17 MeV 4 . Therefore, we can conclude that in presence of a 
strong magnetic field, with trapped neutrinos, the formation of quark matter core 
by deconfinement transition at the central region of a compact neutron star is also 
energetically forbidden. 



4. CONCLUSIONS AND DISCUSSIONS 

We have studied the chemical evolution of quark matter core in presence of a 
strong magnetic field and compared our results with that of field free case. We have 
noticed that in presence of a magnetic field of strength at least greater than B^ e \ 
the critical field strength to populate Landau levels for electrons, quark matter core 
in compact neutron star or proto-neutron star is not possible. The quark matter core 
is energetically unstable with respect to normal hadronic matter in presence of such 
strong magnetic field. On the other hand if a charge neutral SQM in beta equilibrium 
is placed in an external strong magnetic field the system becomes energetically more 
stable compared to field free case. We have also seen that the presence of A-hyperons 
at the core of a neutron star before phase transition to deconfined quark matter 
does not alter the characteristic of newborn quark matter core in presence of strong 
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magnetic field. We have also noticed that our conclusion is insensitive with the 
numerical value of proton fraction. Therefore, the final conclusion is that the quark 
matter core is not at all possible in a compact neutron star or proto-neutron star if 
the magnetic field at the central region exceeds at least the critical value for electron 
to populate Landau levels. In an earlier work we showed that a first order QCD 
phase transition from nuclear matter to quark matter is not possible in presence of 
a strong magnetic field [^jp|,|16|. Even if a metal-insulator type second order phase 
transition is possible, the stability criteria does not allow the formation of quark 
matter core in a compact neutron star or in a proto-neutron star. However, such a 
phase transition is possible in the case of a very old neutron star by matter accretion, 
during slowing down of a pulser which may cause further contraction of neutron star, 
nucleation of SQM by some external sources or by neutrino sparking. Instead of a 
phase transition to quark matter it is possible to have kaon condensed hadronic 
matter at the core region. It would be interesting to see what is the effect of strong 
magnetic field on such condensed state. It can also be a di-quark matter. It would 
also be interesting to study the effect of strong magnetic field on such matter. 
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